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The renormalization of the T-matrix for NN scattering with a contact potential is reexamined 
in a nonperturbative regime through rigorous nonperturbative solutions. Based on the underlying 
theory, it is shown that the ultraviolet divergences in the nonperturbative solutions of the T-matrix 
should be subtracted through "endogenous" counter terms, which in turn leads to a nontrivial 
prescription dependence. Moreover, employing the effective range expansion, the importance of 
imposing physical boundary conditions to remove the nontrivial prescription dependence, especially 
before making physical claims, is discussed and highlighted. As byproducts, some relations between 
the effective range expansion parameters are derived. We also discuss the power counting of the 
couplings for the nucleon-nucleon interactions and other subtle points present within EFT framework 
beyond perturbative treatment. 



INTRODUCTION 

The effective field theory (EFT) method or strategy Q has become a primary tool for studying a variety 
of low energy problems in particle and nuclear physics' important examples include chiral perturbation 
theory (xPT)|2|, heavy quark effective theory (HQET)|3J, and non-relativistic quantum chromodynamics 
(NRQCD)Q. As an EFT parametrizes the short distance physics in a simple way, severe ultraviolet (UV) 
divergences appear. Then one must carefully work out pertinent power counting rules and renormalization 
prescriptions |5|- In nonperturbative regimes (for example, in the application of x?T to low energy nucleon 
systems as is advocated by Weinbergj^), the 'interplay' between power counting schemes and renormalization 
prescriptions becomes quite complicatedM] • To establish a more reasonable and consistent framework, many 
proposals have been put forward 0, H, U E3- creating controversies still to be settled. There were once 
(perhaps are) even doubts about the applicability of the EFT method. 

The main difficulties stem from a distinct feature of the nonperturbative formulation, which invalidates 
the naive use of the perturbative renormalization (through subtraction) programs|lll H^|. In this report, 
we continue our investigations of the renormalization of the EFT for nucleon-nucleon scattering in the 
nonperturbative regimes, which we started in Ref. 11]. We work here with a contact potential that allows 
us to rigorously obtain a closed form of the T-matrix through the use of the Lippmann-Schwinger equation 
(LSE) 13]. In this way, many new features arising in the nonperturbative regime can be explicitly illustrated. 
In our approach, we utilize the underlying theory to understand the renormalization of an EFT. In this work 
we will examine the renormalization prescription dependence of the on-shell T-matrix together with the 
observables or parameters coming from the low energy theorems for nucleon-nucleon scattering. The latter 
are obtained through the effective range expansion This paper is organized as follows: In Sec. II, 
we sketch the nonperturbative parametrization of the T-matrix proposed in Ref.|ll| and its implications 
for nonperturbative renormalization. In Sec. Ill, we employ the algebraic method described in Ref.[13j to 
obtain a rigorous closed form solution of the LSE in the case of contact nucleon-nucleon interactions. Then 
the regularization and renormalization of the T-matrix in the nonperturbative regime are analyzed using the 
closed form solutions at various chiral orders. Both the prescription dependence and its removal from the 
observables or parameters obtained via the effective range expansion (low energy theorems) are investigated 
in Sec. IV. In Sec. V, we study the interplay between the renormalization prescription and the power 
counting schemes for EFT, and the renormalization group evolution in the nonperturbative regime and the 
naturalness of the T-matrix are also explored. Finally, Sec. VI contains our summary. Our main conclusion 
is that one should be aware of the nontrivial renormalization prescription dependence in the nonperturbative 
regime, with emphasis on physical boundaries. 
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A COMPACT PARAMETRIZATION 

Let us start with a standard parametrization for the on-shell partial wave T-matrixQ (we consider the 
diagonal channels for simplicity): 

, 4-7T 1 . . 

Ti-osip =-77 7TT\ ~' 1 

M p cot di (p) — ip 

with M and p being, respectively, the mass and on-shell momentum of a nucleon, and I denoting the angular 
momentum number. Literally, the potential could be systematically constructed or calculated using xPT0 
through counting the powers in terms of p 2 /A 2 or m 2 /A 2 with A being the high scale or upper limit for 
the EFT under consideration (A ~ SOOMeV). Then the off-shell T-matrix for partial wave I could be found 
through the solution of LSE, 

^V l (p',k)Go(k;E+)T l (k,p;E), (2) 

where E + = E + ie with E being the center of mass energy. It should be noted, that Eq- ( EJ) is ill defined, 
as the potential V/ calculated using xPT is usually singular. 

To highlight the nonperturbative features present in this problem, we previously proposed a tentative 
nonpcrturbative parametrization for the T-matrix based on lseHU, 

-Giip,p'\E), (4) 



T^p'-E) Vi(p,p>) 
nl , ^ f0,V l (p' : k)G Q (k;E+)T l (k,p;E) 

Gi(p,p;E) = —^ — — — — — — — , 5 

Vi{p,p')Ti(p,p';E) 

where Qi carries the nonperturbative information of all the quantum processes (all the loop amplitudes in the 
field theoretic terminology) generated by Vi . We have already shown in that the nonperturbative quantity 
Qi could not be renormalized through the introduction of "exogenous" counter terms in the potential 15]. 
That means Q, or equivalently T, is regularization and renormalization (R/R) prescription dependent. There 
can be only one renormalization prescription (contrary to perturbative treatment consistent with physical 
data or boundary conditions. In the following sections we will demonstrate this point with rigorous solutions 
of LSE. 



CONTACT POTENTIAL: RIGOROUS SOLUTIONS 

The specific case of contact potential allows us to transform the integral equation into the algebraic one 
following Ref.^3]. In other words, the nonlocal pion exchange contributions to the potential are neglected. 
However, the main conclusions remain qualitatively valid also when pion exchanges are included. 

Factorized LSE and its algebraic solutions 

Next, to illustrate the role played by the nonperturbative features, we will employ the channel. To 
see how the situation evolves with the chiral orders, we consider the solution of LSE with the following local 
potentials at three different chiral orders (A = 0,2,4) (leading order, next-to-leading order, next-to-next-to- 
leading order): 

A = : Vjf = C ; (6) 
A = 2: V{ 2 f =C + C 2 (p 2 +p' 2 ); (7) 
A = 4 : V^f =Co + C 2 (p 2 +p' 2 ) + C 4P 2 p' 2 + C 4 (p 4 + p' 4 ) . (8) 
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Following Ref.[l3| we "factorize" the potentials into matrices: V — U T XU' with U, and U' being column 
vectors and Aannxn matrix. At next-to- next-to-leading order (A = 4), they are 

Co C2 C4 \ 

C 2 C 4 , U = ( 1, p 2 , p 4 ) , U'=(l, p'\ p' 4 ) ; (9) 
C 4 / 

whereas at next-to-leading order (A = 2) they read, 



(10) 



The off-shell T-matrix factorizes exactly in the same manner: T — U t tII', where r is an n X n matrix. 
Generally, the coupling constants [C„] come from the chiral expansion of an underlying theory (say, QCD) 
for nuclcons and pions in terms of p 2 /A 2 , and hence they scale like C 2 „/Co ~ A~ 2n in the naive power 
counting scheme. 

Using this notation, the LSE can be reduced to the following algebraic equation^], 

t(E + ) = X + Xo1(E + )ot(E + ) (11) 

with o denoting the matrix multiplication. For the 3x3 case, the matrix X and the related parametrizations 
and definitions are listed in Appendix A. With this algebraic parametrization, all the ill defined integrals 
can be isolated and parametrized in any regularization prescription. The solution to this algebraic equation 
is easy to obtain: 

t(B+) = (1-AoI(B+))- 1 oA. (12) 

In a similar fashion, the solution of the T-matrix can be obtained. For the three chiral orders considered so 
far, they read (on-shell), 

A = : -ip, = -J- + Jo + ^-i P ; (13) 

Tosip) Co 47T 
A „ 1 (I-C2J3) 2 r M 

Jos(p) Co + C| J 5 + C 2 (2 - C 2 J 3 )p 2 47r 

1 Ao + A lP 2 + A 2 p 4 M . . 

A = 4: tuF) = d + + d 2P ± + r, 3p6 + Jo + 5F tft (15) 

with the ill defined integrals [J n ] and the coefficients [N n ] and [Z) n ] being defined in Appendix B. The next- 
to-leading order has been considered in Ref.^j|, whereas the A = 4 result has not been given before. We 
should note that Jo always stands "alone" in the real part of the inverse on-shell T-matrix. (A rigorous 
proof of this fact at any order is given in Appendix C.) Here we would like to stress the compact or closed 
form of the expressions for the T-matrix in terms of the couplings and the integrals [J n ]. It is this crucial 
property that distinguishes the nonperturbative solutions from the perturbative ones and complicates the 
renormalization. 

Obviously, the T-matrix becomes more complicated as higher order interactions are included. Neverthe- 
less, when the nonlocal pion exchanges are included, one can naturally anticipate that the nonperturbative 
solutions still take compact or closed forms. Thus we expect, that, our conclusions here will also hold in the 
realistic potentials with nonlocal pion contributions, at least qualitatively. 



Failure of 'exogenous' counter-term renormalization 



Now let us consider the renormalization of the T-matrix. The leading order case is trivial; one can absorb 
the only divergence in Jo into the inverse coupling: 1/Cq, similar to the perturbative cases. 

However, in the presence of higher order interactions, such operation may not work. For example, at 
next-to-leading order, in order to renormalize the on-shell T-matrix in Ea. 114(1 . or to make the fraction 
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Cq+c' 1 Jr +C2(2-C2 J 3 )p 2 ^° nn ^ e ' one should make each of the following compact functions finite at the same 
time: 

(C + C 2 J 5 )/(1-C 2 J 3 ) 2 , C 2 (2-C 2 J 3 )/(1-C 2 J 3 ) 2 , J . (16) 

Now it is clear that the main obstacle for performing "exogenous" subtraction for the T-matrix is the 
compact or closed expressions in terms of [C„] and [J n ]: it is hard to see how to make (Co + Cf Jg)/(1 — 
C2J3) 2 , C%(2 — C 2 ^3)/(l — C 2 J3) 2 and Jo finite simultaneously, since each of them is a compact or closed 
expression expressed given in terms of the two couplings, Co,C 2 , and the three divergent integrals, Jo, J3 
and J5. The situation differs strikingly from the perturbative case where counter terms are introduced order 
by order with the higher order terms discarded, as no compact or 'closed' form of expression is involved. 
Moreover, no matter what was done for (Co + Cf Js)/(1 — C 2 J3) 2 and C 2 (2 — C 2 Jz)/{\ — C 2 J3) 2 , one should 
make sure that Jq stays "separately" finite at the same time. 



Nonperturbative renormalization in EFT 

It is known that an EFT is often established through certain reorganization of parts of a well defined 
underlying theory (UT, at least renormalizable) . Unfortunately, such reorganization usually (1) brings 
about new UV divergences and (2) impedes the exogenous counter terms from working. To see the first 
point, consider the diagrams shown in Figs.l and 2, with the heavy meson exchange diagrams (with g and 
vrih being the coupling constant and the meson mass) underlying the ones with contact interactions. For 
convenience, let us introduce a projection operator P LE to symbolize the influences of this heavy meson: the 
extraction of the EFT vertices or couplings from the UT diagrams. 

At tree level (Fig.l), 

-iC* = V^^=V^{-rp^} = i£ s , (17) 

Ar — m h m h 

no divergence appears. The complication comes at the loop diagram level. For example, in the case of the 
convergent box diagram in Fig. 2, if 7\ E is applied after the loop integration (J ) has been done (the 

correct order), one would get a well defined expansion in terms of -if. When V^e is applied before J , 2 Ai 

h * ' 

(the incorrect order), the divergent bubble diagram results. Thus the new divergences generally arise from 
the incorrect order of computations, as the following commutator does not vanish identically: 



d 4 l 

(2^y 



Oct. =[P™, I 7^4] 7^0- (18) 



Embarrassingly, one has to use EFT either because UT is unavailable or because the calculations in UT are 
tedious. Combining this procedure with nonperturbative context (infinite iteration or resummation) makes 
things even worse: the counter terms could not be implemented exogenously. 

However, from the underlying theory the solution follows immediately: one should devise some procedures 
to effectively "recover" the correct order for P LE and / J§^yr before anything else is done. The clue lies in 
Ea. ((TH|l . Through rearrangement, Ea. ijTB|l is equivalent to the following equation for the integrand of a loop 
diagram (say, the box diagram integrand fbox)> 

d A l f d A l - - f dH 

(2tt) 4 ^ box ~ J (2tt) 4 '^' IjE f box ^~ ^ c • t • f box = J (2ir) 4 ^ bubble ^ c • f box ~ \}-°) 

That means, in order to recover the correct-order results in EFT, we must introduce a counter term: 
Oc .t . fbox- Therefore, the UT scenario provides a natural interpretation for the counter terms, and more 
importantly, a rationality for the subtraction at the level of the loop integral without any reference to La- 
grangian. That is, the counter terms must be endogenous: the divergent integrals must be subtracted before 
the nonperturbative reorganization^]] ■ Finally, the subtracted integrals (finite) appear in the compact non- 
perturbative expressions which are no longer compatible with exogenous counter terms. In this logic, the 
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formal consistency issue0 of the Weinberg power counting simply dissolves: in the nonperturbative regime 
of EFT, there is no point in searching for exogenous counter terms and their counting rules. Of course there 
might be other approaches that directly renormalize the integrals without the explicit use of counter terms. 

Thus the nonperturbative renormalization must be implemented either through endogenous counter terms 
or other means that effectively "subtract" the EFT integrals (see similar prescriptions in Ref.^3)- For the T- 
matrix considered in this paper, this procedure would be formulated as a simple replacement of the divergent 
[ J n ] with the subtracted [J n ] , which are finite constants (prescription dependent and hence arbitrary) , in the 
compact expressions: 

A = : f^ 1 (p) = ±- + J + M- ip - (20) 

A = 2 : f °° {p) = c + c^c 2 f-c 2 j^ + J ° + ^ 

A=4:T - (p) = ^www +Jo+ ^- (22) 

Here we wish to note that from a UT perspective, both [C n ] and [J n ] come from the projection acting on 
the convergent vertices in UT. (Those that are divergent in UT will be renormalized before applying the 
projection and do not directly contribute to EFT renormalization due to scale hierarchy, we will return later 
to this point in Sec. VI.) In this sense, [J n ] are also fundamental parameters in EFT, so the nonperturbative 
T-matrix is parametrized by both [C n ] and [</„]. To illustrate this point, let us apply the projection on the 
box diagram after the loop integration is carried out. After some calculations we get 

^ifilleadingterm = -i^I^M, ™ h ,p) = iC*[J^(M,m h ) + ^ip], (23) 

where the definite parameter Jg UT) (M, m/,) (see Appendix B), in a place of the divergent integral Jo, can be 
extracted in the following way 



\1f le! tree) (. \ lL tree) 



jr\M,m h ) = - ' „ U) b °:: \ P = = -v^ i P=0 . (24) 



We should note that here we used the mesonic interaction for illustration. Of course, the true contact 
nucleon interactions should be computed from QCD. But the mechanism explained above still holds in 
general. 

Finally, we note that the renormalized T-matrix suffers from severe prescription dependence in the non- 
perturbative regime, which is incompatible with the exogenous counter terms. That means, given specific 
couplings, only one prescription could yield the physical T-matrix; others have to be dropped even though 
they are finite. So the final resolution boils down to the flexible regularization methods that could facilitate 
convenient access to physical predictions |ll|. as was already noted in other nonperturbative contexts|l8|. 
This argument leads us to the following strategy: One first parametrizes the ill defined integrals in terms 
of ambiguous constants and then imposes physical boundary conditions. Similar strategy also based on the 
underlying theory, has already been described in Ref.|l9| for renormalizing any EFT. 

LOW ENERGY THEOREMS (LET) AND PRESCRIPTION DEPENDENCE 

Effective range expansion 

Now let us consider effective range expansion (ERE) defined as follows, 

{4 1 11 x 

- jjJ^GO \ =PCotS(p) =-- + -r e p 2 + VkP 2k , (25) 
' k=2 



G 



with the parameters a and r e being the scattering length and the effective range, which (including [vk]) 
could be extracted from the scattering data. In this sense, we could impose their values as the boundary 
conditions for the T-matrix. Performing the expansion for the T-matrix obtained above, we get: 

A = 0:pcot%) =-_{C - 1 + Jo}; (26) 
A = 2 : p cot S(p) =~\ voK 1 + Jo ~ *<M~ V + £ Mi *o k_1 (V)* \ , 



k=2 

\2 <F — r< i r<2 



D =(l-C 2 J 3 y, 5 Q = C + CiJ 5 , S 1 = C 2 (2 - C 2 J 3 ): (27) 

4-7T - - - — - — — 

A = 4 : pcot J(p) = - — {NoD^ 1 + (NtDo - N^D^p 2 

+ [N 2 D 2 - JViDiA) + N (Dl - D D 2 )}D^p 4 

+ [^0(2/3x732/30 - D 3 D 2 - Df) + NMD 2 - D D 2 ) - N^D^D^p 6 

+ •••}. (28) 

The scattering length, effective range, and the v k can be read from the results above. For the three orders 
considered so far, we have 

A = 0: a- 1 = ^(C^ + Jo), r e = 0, v k =0,k>2; (29) 
A = 2: a-^gfa^ + Jb), r e = ^u ~8^ 2 , v k = ^v 5^-5o)~ k ~ 1 , k > 2; (30) 
A = 4: a- l = ^{N D^ + J ), r e = ^(NoD, ~ N 1 D )D 2 , 
v 2 = —[N (D D 2 - D 2 ) + N\DiDq - N 2 D 2 ]D^, 

A'jr — — — — 

v 3 = —[N (D 3 D 2 - 2D Y b 2 bv + Df) - JViA>(5? ~ D D 2 ) + N 2 DiD 2 ]Dq 4 , (31) 

Note again that, at each order, Jo only enters the expression for the scattering length but is "decoupled" 
with all the other ERE parameters. The reason is clear: Jo stands alone in T^ 1 . As we make clear below, 
this point has very important implications. 

Keeping in mind that the parameters [J n ] are in principle independent of each other, two distinct ap- 
proaches can be adopted in order to impose physical or reasonable boundary conditions: (1) taking the 
couplings and the prescription parameters as the fundamental variables, and the ERE parameters as the 
functions of these variables; or (2) conversely, taking some of the ERE parameters (which should be physical) 
as fundamental and the others as the functions of them. For convenience, we could also parametrize [J n ] in 
terms of a dimensional scale Jl and dimensionless numbers [<?...]: 

Jo = qoMJl, J 3 = q 3 MJl 3 , J 5 = q 5 MJl 5 , J 7 = q 7 MJl 7 , J 9 = qgMJi 9 . (32) 

The appearance of M is easy to see from Appendix A. Generally, the magnitude of Jl could vary from a value 
on the EFT expansion scale, A, to the value of the pion decay constant (much smaller than M): Jl e (f v ,k). 
One can also alter the integrals by letting the dimensionless numbers [q n ] to vary. Thus, a nonperturbative 
rcnormalization prescription is parametrized by [q n ;Jl]. However, the magnitude of [J n ], which also comes 
from the low energy projection in UT, should not be larger than the naive powers of the chiral symmetry 
breaking scale A x sb — M. In other words, we can safely assume that: | J„| < M n+1 ,n ^ 0, | Jo| < M 2 . 

Having made these preparations, we can start to examine the low energy expansions listed above order 
by order in chiral expansion. The leading order case (A = 0) is trivial: we have only one condition, i.e., 
imposing that ^ + Jo = a, with a being experimentally measured, is enough since r e = v k = 0, k > 2, which 
is obviously bad theoretical prediction, although it is not prescription dependent. Thus, the situation at this 
order is physically uninteresting, and most importantly, the distinctive nonperturbative features we wish to 
expose are not obvious here. Therefore we examine in detail the higher order cases. 
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LET at next-to-leading order: A = 2 

Let us start with the next-to-leading order: A = 2. As was mentioned before, we shall discuss the problem 
from two perspectives, respectively. 

First point of view 

For convenience, let us list the explicit expressions of a,r e , etc. in terms of [C n ] and [J n ], which read, 

_ M Cp + Clh 

a ~ 4tt (1 - C 2 J 3 ) 2 + Jo(C + Cf J 5 ) ' 

= 8^ (2C 2 ~Clh)(\-C 2 ,hf 
Te M (Co + C 2 2 J 5 ) 2 

fe "M () (c + cp 5 ) k+1 '"" 2 - (33) 

Before imposing reasonable boundary conditions, we could not make any physical predictions as Jo, J3 
and J5 each independently varies. The independent variations of Jo, J3 and J5 could not be easily absorbed 
into the couplings due to the reasons explained in Sec. III. Thus, unlike the leading order case, we need 
(more) exogenous constraints to fix the values of Jo, J3 and J5. 

From Ea. (|33[l . J3 and J5 could be solved in terms of r e , v 2 and Co, C2, the solution might be unique after 
accounting for a reasonable size. After the insertion of the obtained numbers back into the formula for a, Jo 
could be expressed in terms of the physical value of the scattering length. In this sense, imposing two addi- 
tional boundary conditions could fix the prescription or make the next-to- leading order result unambiguous. 
Now the [vk] with k > 3 are taken to be theoretical predictions, which now shall be better than the leading 
order ones, as we have more degrees of freedom to work with: J3 and J5, which come together with the new 
interactions. The fact, that some of the predictions are still poor can be attributed to the inadequacy of 
the next-to-leading order potential: even higher order terms should be put in and accordingly more phys- 
ical boundary conditions are needed. Thus, in spite of the fact, that the procedure for fixing prescription 
becomes more nontrivial and laborious, the predictions for the ERE parameters improve when the higher 
order interactions are included (because of the "freedom" brought by the augmented interactions). We must 
repeat here that the predictions are made using the prescription that is most compatible with the physical 
boundary conditions. 

Employing 1|32[1 . we could also write the above equations as 

_ M Co + Clq 5 Mfi 5 

° ~ 4tt (1 - C 2 q z M^f + q MJi(C + C$qiMJP) ' 
_ 8^ (2C 2 - CjggMgXl - C 2 g i Mp i ) 2 

Te ~ M {C Q + cq q5 Mji 5 ) 2 

(C + Clq 5 MF) k+1 

Note that even though the Jl scale dependence can be removed (fixed), the prescription dependence remains 
in terms of the dimensionless parameters [q n ] that are independent of each other, a subtle point that often 
seems to be overlooked. For example, in a cutoff scheme, the renormalization is usually performed in such a 
way that the cutoff dependence is removed by letting the couplings to develop a certain cutoff dependence. 
However, remains residual prescription dependence remains through the cutoff independent but prescription 
dependent numbers [qn]- Without fully appreciating this point, any fitting procedure that uses tuning of the 
couplings or even tuning of the cutoff only amounts to fitting along a special orbit in the space [q n ] , not in 
the whole space. The result thus obtained is still prescription dependent. 
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Second perspective 



Taking a, r e and {^2} as the elementary parameters, we can express all the higher order constants [vk, k > 3] 

as 

v k = ^^(-2 U2 /r e ) fe ,fc>3. (35) 

Note that from this perspective the prescription-dependent parameter Jo seems to be an independent constant 
in addition to the three elementary parameters. At first sight, this ambiguity calls for one more condition: 
the value for v 3 . But from the discussion above, we know that this seemingly independent parameter is in 
fact determined together with J3 and J5 by the equations l|33|) ■ Of course the nature of the problem remains 
the same even when it is taken as independent. The most striking point here is again that the prescription is 
"removed" through fixing, i.e., through boundary conditions, which is a nontrivial procedure as articulated 
above. 

One could also put Ea. (|35Jl into the form that contains no explicit prescription dependence: 

v k =v 3 {-2v 2 /r e ) k ~ 3 ,k>4. (36) 

This relation should hold for any problems (in certain atomic or molecular contexts) with contact potential 
such as V(x) ~ C* (0)( 5( 3 )(x) + C* (2) V 2 ^ 3 )(x). 

So, we may conclude that, no matter what point of view one adopts, the key point is that the compact 
nonperturbative formulations make the prescription dependence and its removal a very nontrivial problem 
or procedure. However, the predictions also improve with the use of these formulations, in spite of their 
technical complexity. This nontrivial procedure will get more involved as more higher order corrections to 
the potential are included. To verify this, let us turn to the next-to-next-to-leading order. 



LET at next-to-next-to-leading order: A = 4 



Again we begin with the first perspective. 



First perspective 



Now we have four couplings (Co, C 2 , C4, C4) and five prescription dependent parameters (Jo, J3, J5, J7, Jg) 
in eight compact expressions: No, N\, N 2 , D , D\, D 2 , D 3 and Jo, which stays alone. From Eq.J2U it is clear 
that we need at least five conditions to fix Jo, J3, J5, J7 and Jg, say a,r e ,v 2 ,v 3 and U4. But the compact 
expressions such as that for the scattering length, a = |^ — ®° D with Dq and No given in Appendix B, 

become more involved. This means that the boundary conditions might be more stringent for Jo, J3, J5, J7 
and Jg and the analytical work more difficult. In the meantime, the predictions for Vk at this order should 
be better than the leading and next-to-leading orders, as we have more parameters. 

Here some remarks are in order. At next-to-leading order, we ignored the possible multiple solutions for 
the fixing procedure. Here, with more compact expressions being involved, we should be more careful about 
this multiplicity of solutions. To this end, we note that the multiplicity could be effectively reduced with 
the limitations on the reasonable magnitudes of [J n ] together with the experimental values of the higher 
ERE parameters (say, Vk,k > 3). However, no matter how the multiplicity is removed, the solution is still 
an approximate (though nonperturbative) one: The equations (I29I3UI3I|I are obtained from a truncated 
potential and could not be exact ones. Then, the theoretical predictions based on such equations will be less 
credible, especially for the ERE parameters that dominate higher and higher energy regions. In other words, 
the boundary conditions should be given by a procedure similar to fitting the shape of the phase shift within 
the corresponding ranges at each chiral order. This is actually what most authors have done, though the 
rcgularization schemes used vary significantly. Of course our remarks in Sec. IV.B.l concerning the residual 
prescription dependence still apply for all the higher order calculations. 

Mathematically, the multiplicity of solutions might be generic for nonperturbative renormalization because 
of the compact expressions involved. Therefore the limit cycles encountered in the Schrodinger approach 
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of renormalizing singular potentials |20| might just be examples of such multiplicity in certain regularization 
schemes. 

Second perspective 

To discuss the problem from the second perspective, we need to express everything in terms of the first 
five ERE parameters. Then the arguments go as in the preceding subsection. We shall not, however, repeat 
such complicated technical details here. One could also find the relations like Eq.(35) or (36), which hold 
true independently of the prescriptions, by repeatedly using the following recursive relations: 

Dk As _i D 2 Di D 3 D 2 Di 

v n = - > -=-?;„_fc,n > 5; v 3 = -=-a - T ~r e - -=-v 2 ; v± = --p-r e - -e-v 2 - -=-v 3 ; (37) 
Uq JJq LJ D Uq ij l> 

_i M - M M 

a = - J , r e = —r e , v n =—v n ,n>2, (38) 

where the coefficients [j^] could be solved in terms of M, a, r e ,v 2 , v 3 and V4,. 

Lessons from nonperturbative solutions 

Now it is clear that things get more complicated as more higher order terms are included in the potential. 
Given this, we should not take the lower order results too seriously. For instance, the low energy theorems at 
leading order are too simple to be true in practice: r e = Vk = 0, k > 2. This implies the necessity to include 
higher order terms, which, however, will bring us both favorable and unfavorable consequences. On one hand, 
more severe prescription dependence will show up and make this analysis more difficult. On the other hand, 
more prescription ambiguities also provide us with more chances to access the measured values of the ERE 
parameters. Although our calculations were done for the case of contact potentials, the core feature of our 
analysis-more ambiguities or more divergences at higher orders-holds true also for realistic potentials. At this 
stage, we shall mention that the freedoms in the prescription are in fact limited: [J n ] must satisfy certain 
requirements as presented in the discussion following Ea. (|32[) . Moreover, the coupling constants should 
generally follow certain rules of EFT power counting. Then, after putting all these theoretical aspects into 
consideration, the EFT predictions must lie in certain region of the 'space' of observables. 

Now we provide another way to see the virtue of the fitting procedure. Let us examine the variation of 
the functional form of the scattering length a in terms of the couplings and [J n ] for different chiral orders: 

M Co 



A = 



A = 4 



47rl + C Jo' 

M Cp + Clh 

4tt (1 - C 2 J 3 ) 2 + Jo(Co+C 2 2 J 5 )' 

M Do(Co, C2, C4, C4; J5, J7, Jg) 



47r N (C , C 2 ,Ci, C4; J 3 , J 5 , J7, J 9 ) + J a D (Co, C 2 , C 4 , C 4 ; J 5 , J7, J 9 ) 

It is obvious that the theoretical form of the scattering length varies with the chiral order quite significantly! 
So the scattering length calculated at lower orders should not be directly identified with the experimental 
value in order to accommodate the higher order terms. Thus a more reasonable way for fixing the renormal- 
ization prescription will be the one that avoids the direct identification of physical parameters in order to 
accommodate higher order contributions in a consistent way. To this end, again a procedure like fitting the 
empirical curve over appropriate low energy regions might be more plausible. 

POWER COUNTING AND RENORMALIZATION IN THE NONPERTURBATIVE REGIME 



In all the discussions above, we have left out the power counting of the couplings. Since they constitute 
the basis for the EFT methods, it is necessary to see what the nonperturbative renormalization procedure 
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described above means for the power counting rules. In fact, as was stressed in Ref.^jJ, the parametrization 
in the nonperturbative regime given in Eq.Q implies that, in order for a power counting scheme for cou- 
plings to be meaningful, the corresponding prescription for the constants [J n ] must be appropriately chosen. 
Otherwise, one could not obtain the physical T-matrix. 

From the standpoint of UT, both [C n ] and [J n ] come from the well defined low energy projection (P LE ) 
applied to UT amplitudes. So both [C n ] and [J n ] serve as the elementary parameters for parametrizing the 
T-matrix for the low energy nucleon-nucleon scattering. In the EFT treatment without knowledge of the 
details from UT, we are forced to employ [C n ] as the elementary couplings according to certain counting 
rules, while the constants [J n ] appear as the divergent pieces in the EFT loops constructed with the use of 
[C n ]- Thus it is the EFT treatment that makes [C n ] and [J n ] look disparate. In UT they are organized and 
derived together according to more elementary rules. Therefore, changing any of them (each single parameter 
in [C n ] U[J n ]) alone would alter the physical behavior of T. Thus they must be considered together. 

One could also understand it from the Wilsonian definition of EFT through successive decimation of the 
higher scales, where different EFT expansion point would lead to both different couplings and different [Jn], 
as long as the expansions are compatible with the chiral power counting. 

To be specific, the variations of [C n ] and [J n ] (from now on Jo is excluded from [J...] for the reasons to be 
given below) must not alter the functional form (shape) of the T-matrix: 

Re f i _ T \ = Emc:..];[j'j) P 21 = Emc...];[j..])p 21 = Eivfoy* m] 

\T os (p) °J EA(P'.J;R'J)P 23 V /),: Y'... : ■/...!/;-'■' v/;;v v * 1 1 

Here we use the superscript 'phys' to indicate that the parameters in the last fraction are physically de- 
termined, for example, from a genuine UT. To see why J is excluded from [J n ], consider the physical 
parametrization of T-matrix (independent of the variations of [C...] and [Jn]), which has the following form 

- ^ 1 1 ' Mj + ^-ip, (40) 



where 7 must be a physical scale, just like the nucleon mass M and the on-shell momentum p. Now it 
is clear that Jo alone corresponds to the physical parameter M7, which should therefore be independent 
of prescriptions. If it were not the case, or if Jo could vary with prescriptions, we would have to alter 
Nq and Dq to compensate for such variation in Jo. Now, to keep the proportionality between Nq,Dq and 
NiP 2t , Djp 2 ^ invariant, all the rest of [N,.., D...] must be accordingly altered, which in turn leads to an overall 
factor for A^'V'VCE p 2 ?). Then, the functional dependence of the T-matrix upon p would be 
altered, since its imaginary part, ¥-ip, remains intact. Hence Jo must stay independent of prescriptions, i.e., 
physical. One could also verify this by examining the consequences on the ERE parameters. 

As p is arbitrary in the supposed range, the 'invariance' discussed above leads to the following nontrivial 
equations for [C n ] and [J n ] with the crucial presence of the physical parameters [N- ?h7 "^] and [Z)j phrs ' 1 ] for the 
on-shell T-matrix, 

Ni(P...]\ [J J) = N ( r\ ^([C.J; [J J) = DW, (41) 

These equations have dual implications: they could be used either (1) to fix the prescription ([J...]) in terms of 
the couplings ([C...]) and the physical parameters ([_/V( phys ); -D^"']) or conversely (2) to examine the influence 
of prescription upon the couplings with the help of the physical parameters. The first use just parallels what 
we have done in Sec. IV. B and C. 



Interplay between Power counting and prescription: next-to-leading order 

Let us illustrate the interplay between power counting and prescription at next-to-leading order; that is, 
we try to solve the following equations for couplings: 

= «0 = — 7— T, ~7Z 7T^="2 = — 7— T . (42) 



(l-C 2 J 3 ) 2 U_ A^)' (l-C 2 J 3 ) 2 '- N ^> 
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The solutions are easy to find as, 

C^=J3 1 {l±(l + a 2 Ja)-^}, (43) 

^ ±) = T ^-|{l±(lWa)-} 2 . (44) 
i + a 2 j 3 J 3 

Taking into account the natural boundary condition for C 2 : C 2 \j^q =>• & 2 /2, wc are left with unique 
solution: (and Cq - '). Thus assigning a power counting to Co and C 2 means assigning the sophisticated 

scaling for J 3 and J5. Conversely, one can come up with an alternative interpretation: The power counting 
for the couplings could only be preserved in some particular prescription in order to obtain the expected 
physical behavior from the T-matrix. Note that here we have deliberately not mentioned Jo, it will be 
exclusively discussed below. Equations l|41|) or (|42[1 now formalize our discussions concerning the interplay 
between power counting and prescription. 

More interestingly, these equations have a further utility: they could be used to describe the evolution of 
the couplings in terms of a sliding scale (/i) in [J n }{— [q n Mfj, n ,n ^ 0]). Since the exogenous counter terms 
are incompatible with the closed form of the T-matrix, the conventional route to the evolution described by 
renormalization group equation does not exist. But we could take the evolution implied by Eqs. (|41|1 or (|42|) 
as a nonperturbative "renormalization group" evolution. We discuss this point in the next subsection. 



Nonperturbative 'renormalization group' (RG) evolution 

To proceed, let us choose the prescription with [J n = q n M/i n ] to examine the evolution of the couplings 
enforced by l|42|l . Let us assume that there exist enough boundary conditions to obtain the 'physical' solutions 
for the couplings from the equations in l|41|) : 

d = Fi([N<*»\D<*r\M\\ [g...];/*),Vi. (45) 

With such nonperturbative solutions, the complete evolution of the couplings are determined and both the IR 
and the UV fixed points can be identified. For example, at next-to-leading order, we have from Eas. (|44l43|l . 

C (ao, a 2 ,M, q 3 , q 5 ;^) = "° - -£-{1 - (1 + q 3 a 2 Mfx 3 )-i} 2 , (46) 

1 + q 3 a 2 M(J, d q^M/j, 

C 2 (a ,a 2 ,M,q 3 ,q 5 ;n) = {q 3 M^)~ x {l - (1 + q 3 a 2 Al n 3 )~i}. (47) 
It is easy to see that they have both IR and UV fixed points: 

IR fixed pointy =*> 0) : Ctf R) = a , C ( 2 IR) = a 2 /2; (48) 
UV fixed pointy =*> 00) : C^ UV) = cf V) = 0. (49) 

Note that the prescription dependence is obvious in Ea.(|45|) with the presence of [q...], but the UV and 
IR fixed points are prescription independent. While the IR fixed points are realistic as the couplings were 
defined in the low energy limit, the UV fixed points seem not to be realistic. But such UV behavior of the 
EFT couplings is compatible with the fact that the EFT couplings would be dominated by the UT couplings 
at high energy, and therefore 'vanish'. Of course we should bear in mind that, what we obtained are only 
approximate answers, though nonperturbative. 

Note that the Eas. (|41l45(l contain the full dependence upon the prescription parameters. So one could 
also derive the equations a la Stiickelberg and Petermann[2l| that describe the laws for transitions from one 
prescription to another which are not related by running the renormalization scale: 



{Ni,Dj} =0,Wi,j. (50) 



In terms of [q...; fj] they become 



d[J. 

^{JV i ,D,}=0,Vi ! i. (51) 
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In the foregoing discussions the physical requirements arc imposed on the functional shape of the T-matrix. 
Alternatively, we could also employ the physically determined ERE parameters (scattering length, effective 
range, etc) instead of [N^**'; D^*'] to solve the couplings in terms of [J n j. In principle, the two approaches 
should lead to the same evolution behavior, but the ERE approach is more involved than the shape approach 
as is clear from the comparison between the Eqs. <|3(JI31|) and Eas. (|41|l . 



Determinant for the natural or unnatural scattering length: Jo 

Now we discuss the determinant (s) of the size of physical parameters [N- phy ^] and [T)j phys * 1 ] or a, r e , Vk, Vfc > 
2. As argued above, a complete parametrization of the T-matrix is given by [C...] supplemented with [J...], 
then [C..J and [J...] together determine whether the physical parameters are of natural size or not. Wc 
could have four rough scenarios, listed in Table [I] where by a natural C n we mean that the scale A in its 
parametrization C n ~ 1/(MA" +1 ) is of the size of the expansion scale (unnatural if A ~ p, m^), while for 
[J... = q.„Mfi'"] the situation is reversed: the natural size of fi should be ~ p, m^. A natural T-matrix is 
parametrized by [A^ phrs ' ) , D^^] (or for a,r e ,Vk,^k > 2) such that the dimensional parameters are of the 
same magnitudes as the natural couplings. 

Examining the concrete expressions of the T-matrix, we find that whether the T-matrix is natural or not 
is determined by both the sizes of the couplings and the magnitudes of the dimensionless combinations like 
Iln.m^n 1 '^! ( dim [EL,m 1 Jm] = 0). Now suppose we have natural couplings, i.e., C n ~ l/(MA n+1 ), Vn. 
If [J...] are also natural, we should have | J\ m C„J m \ <C 1 for all the dimensionless combinations. Then, 
given our experience at next-to- leading and next-to- next-to-leading orders, we can anticipate that: 

N^~l, N t ™~-L; D ™~ 1 J—(~C j ),ViJ. (52) 

In such case we shall obtain a natural T-matrix, or ERE parameters (a,r e ,ffc,Vfc > 2), of natural sizes. If 
[J...] are unnatural, then in general, we could have | m C^ 1 J m \ ~ 1 for the dimensionless combinations. 
Therefore, we have, 



A^'»l(«l), iV<^> » («)-!-; ^>>( < <)_i_ ( ^Q),V i ,j. (53) 



In this case, we obtain an unnatural T-matrix, or unnatural ERE parameters (o, r e , Vk, Vfc > 2) with natural 
couplings. For example, at next-to-leading order, we have, 

natural [J..] : \C 2 J 3 \ < 1, |Cf <V J 5 | < 1, 

^2 



natural T : ~ ■ t I 7i - r^or^ (54) 



(1 - C 2 J 3 f l 

Co + ClJ h + C 2 (2- C 2 J 3 )p 2 ~ Co + 2C 2P 
unnatural [J..] : \C 2 J 3 \ ~ 1, ICfC^ 1 J 5 \ ~ 1 



(1 — C 2 J 3 ) Ci 
=> unnatural T : = ______ (55) 

where each of can be either pretty small or pretty large and therefore the T-matrix could not be a natural 
one. 

Now let us consider Jo and the scattering length in particular. As argued above, Jo should be viewed as 
an independent physical parameter, not as a common prescription parameter. From the parametrization of 
T and the formulae in the preceding sections, Jo will only contribute to the scattering length: 

-i M- N ([C..,JJ) 
a — — Wo 



M \ D ([C...,JJ) 
8^ (N ([C...,J.])D 1 ([C...,JJ) ^([C.., J...]) 



M{ f)2([C... ! J...]) D ([C..,JJ) 
da^ 1 An dr e dvk 



9 J M ' d Jo dJ 



0,Vfc>2. (56) 
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Now we can see that, even when both [C..J and [J...] are of natural sizes, the scattering length could be 
unnaturally large once Jo is unnatural (~ MA), 

naturalJ (- Mfx) : a" 1 ~ -0(A) + 0(/i) - -0(A); (57) 
unnaturalJ (- MA) : cU 1 ~ -0(A) + 0(A) 0(/i). (58) 

That is, in the 1 S'o channel, there theoretically exists such a scenario that the scattering length could be 
unnaturally large while all the rest ERE parameters arc naturally sized. Then the first situation in Table Q] 
should be amended as follows: even when all the couplings and all the rest [J.J are natural, we would get 
an unnatural scattering length as long as Jo is unnatural. 

For an unnatural power counting of the couplings, the discussion would be more difficult and we refrain 
from exploring such situations here. As we have shown, that both the natural and the unnatural physical 
parameters could be explained with the natural couplings (provided the nontrivial nonperturbative prescrip- 
tion dependence is fully explored), we feel that it is more reasonable to work with natural or conventional 
power counting of EFT couplings. 

DISCUSSIONS AND SUMMARY 

Now it is time for us to address some theoretical aspects that have been omitted or not fully discussed so far. 
Let us start with the relation between the UT renormalization and the EFT renormalization that is involved 
in Sec. III.C. Generally, in an EFT one deals with the new divergences in the diagrams that are induced 
by some low energy expansion (with the "wrong" order of operations, C.f. Sec. III.C.) or similar operation 
in UT, whereas the diagrams that need renormalization in UT are usually hidden in the EFT couplings. 
Note, that the diagrams that renormalize UT dominate the quantum fluctuations at short distances, while 
the ones that are divergent in EFT dominate those at long distances. So the two renormalization do not 
interfere with (affect) each other due to the large scale hierarchy between UT and EFT, i.e., they work at 
two widely separated scales. Thus, the renormalization in UT does not affect the renormalization in EFT. 
This supplements our remarks after Ea. (|22[) in Sec. III.C. 

Next, let us address the effect of the potential truncation on the nonperturbative renormalization group 
evolution. At any fixed chiral order, the nonperturbative evolution behaviors of the coupling of the highest 
chiral dimension should be less trustworthy. This is because once the next order interactions are included, 
the coefficients for the term with the highest power of p would suffer the largest changes in the functional 
forms; the coefficients for the lower power terms receive smaller changes from the new couplings. That means 
that due to the truncation of the potential, the nonperturbative evolution behaviors of the EFT couplings 
with lower chiral dimensions should be more trustworthy than those with higher chiral dimensions. One 
could see this point by noting how the forms of [JVj([C...]; [J..]), Di([C...]; [J..])] (as functions of the couplings 
[C...]) change with the inclusion of higher order interactions. 

In Sec. V.C., we have shown that a natural (or conventional) chiral power counting of the EFT couplings 
does allow the T-matrix to have unnatural parameters, or unnatural scattering length, etc. In particular, 
there is a possibility that only the scattering length is unnatural while the rest of the parameters are natural. 
This seems to be just the realistic situation with the 1 S'o channel nucleon-nucleon scattering at low energy. 
This scenario is clearly different from the one discussed in the literature where unusual power counting of the 
couplings was employed 22] . Here the key role is played by the nonperturbative renormalization prescription. 

Although our conclusions or remarks have been reached with contact interactions, we feel that the con- 
clusions or scenarios depicted here should remain qualitatively true even in a realistic situation because the 
crucial features of the nonperturbative renormalization remain unchanged: (1) More ill defined pieces in the 
loop integrals appear at higher chiral orders; (2) The nonperturbative solution of the T-matrix takes a closed 
form that can only be renormalized via endogenous counter terms. Alternatively, one could also take the 
rational function form as a Pade approximant to the realistic T-matrix. 

Now let us comment on the literature. In Ref.[2j|, a subtraction similar to the endogenous one described in 
the present paper is employed: the counter term is introduced before the T-matrix is calculated, a procedure 
that parallels the loop integrations. However, it is not clear if the subtraction described in some pap ers is 
equivalent to the endogenous one or not. For example, the subtraction procedure described in Ref.[24( does 
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TABLE I: Determinant for naturalness/unnaturalness of T-matrix 





Natural [C...] 


Unnatural [C...] 


Natural [J...] 


Natural T-matrix 


Unnatural T-matrix 


Unnatural [J...] 


Unnatural T-matrix 


Natural/Unnatural T ? 



not appear to be an endogenous one. Thus it may be flawed, as was already noted in Rcf. 25] . In Rcf. 25], 
the whole investigation is made in the nonperturbative formulation (compact) of the T-matrix, a positive 
aspect of this study. However, a special regularization (cutoff regularization) exclusively used in Ref.|25l] 
unfortunately makes their analysis inevitably prescription dependent. In a contrast, the strategy employed 
in Ref. 9] for parametrizing and fixing the nonperturbative renormalization prescription dependence is closer 
to the one used in the present paper. The importance of boundary conditions has already been stressed in 
Ref.|ld|. where the physical observables, such as phase shifts, were parametrized without involving explicit 
divergences. 

Obviously, we just explored some convenient scenarios of the nonperturbative solutions. Our arguments 
have been unable to exclude many other possible scenarios. The only point in favor of the scenarios discussed 
in this paper is that they are relatively simple, whereas the rest possibilities seem rather sophisticated, and 
often use fine tuning or similar arguments. 

In our opinion, a better way to work with the renormalization in the nonperturbative regime is to appreciate 
the presence of a well-defined theory underlying an EFT, as illustrated in this paper. In this sense, the 
renormalization of singular potentials in quantum mechanics |26f. or equivalently the self-adjoint extension 
of singular operators in Hilbert space, should also be embedded in the underlying theory background. This 
is plausible since quantum mechanics IS an effective theory of quantum field theory. 

In summary, we reconsidered the renormalization of the EFT for nucleon-nucleon scattering in the nonper- 
turbative regime using contact potentials that facilitate rigorous solutions of LSE. Detailed analysis reveals 
that the T-matrix in the nonperturbative regime should be renormalized through the endogenous counter 
terms whose net effects are to remove the divergences in the loop integrals, or through means that could 
yield the same results. The rationality for the subtractions at loop integral level could be naturally explained 
from the underlying theory, with the UV divergences being shown to come from the "incorrect" order of 
operations in the construction of EFT. Then, using the effective range expansion, we demonstrated that the 
nontrivial renormalization prescription dependence in the nonperturbative regime must be "removed" by 
imposing appropriate boundary conditions. We also argued that when imposing boundary conditions, the 
full " space" for renormalization prescriptions should be explored in order to be able to remove any residual 
prescription dependence. It is also important to impose the boundary conditions in such a way that higher 
order terms in the potential could be consistently incorporated. Finally, the nontrivial relation between the 
power counting of the couplings and the renormalization prescription was highlighted in the nonperturbative 
regime. As byproducts, (1) the nonperturbative 'renormalization group' evolution was described; (2) the 
naturalness of the scattering length, etc. were shown to be compatible with the natural or conventional 
power counting of the couplings because of the nontrivial prescription dependence. That is, the nontrivial 
prescription dependence becomes a virtue in such a case. Obviously, much work remains to be done. 
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Appendix A 



/ f d 3 k : r d 3 k fc 2 



Z(25+) 



T2¥F ' J T3FF ■ 

f ri 3 fc fc 4 



f dVt fc 

J J2^Y E+- 



f d 3 k k l \ 
f d 3 k k e 



\ 



r d"fc fc* r o fc 

The entries of this matrix can be parametrized as follows (p = \JME) : 

.Mp 



E+-l§ J 



I 



I 



d 3 k 


1 




(2tt)3 


£+ - 


fc 2 

M 


d 3 k 


e 




(2tt) 3 


E+ - 


fc 2 
M 


d 3 k 


fc 4 




(2tt) 3 


E+ - 


fc 2 
M 


d 3 k 


k e 




(2tt) 3 


E+ - 


fc 2 
M 


rf 3 /c 


fc 8 





— — Jo - 

J 3 + %> 2 ; 



47T 



(27T) 3 E+-% 



■k + -hp 2 + hp A \ 



Ji + hp 2 + J 3 / + /o/; 

Jo + J7P 2 + J 5 P 4 + J3P 6 + hp 8 - 



(59) 



(60) 
(61) 
(62) 
(63) 
(64) 



Here { J„} with n — 0,3, 5, 7, 9 are regularization and renormalization prescription dependent constants. 



Appendix B 



N = 

m = 

N 2 = 
D Q = 
Dx = 
D 2 = 
D 3 = 



J< UT) (M,m,) 



(1 — C 2 J3 — C4 J5) 2 — C0C4J3 

— C^Cf jf Jg + 2C4C4J3J5J7; 



2C4J3 

-2. 
3 ' 



C4J5 + 2C4C4J5 — C4C4 J 5 — 2C4C4J3J7 
C*4 J3 + 2C2C4 j| + 2C4C4 J3 J5 + 2C| J3 J5 — C4C4 J3 j| + C4CI J3 J7; 

Co + C| J5 + Cf J9 — C0C4J5 + C4C4 J 2 + 2C2C4J7 — Cf C4J5 Jg; 

4 J7 + 2C4C4JJ — C4C4J5J7 + 
2C2C4J3 — 2C4C4J5 — C 4 J$ + C4C4 J5 — C4C4 J3 J7; 

(y + 3 - 2a;) 2 + 8(a: - 1) 



2C2 — C 2 J3 + C0C4 J3 



2C4 + C4 

y 4 



-C\J%. 



mi 



(4n) 2 M 2 



f dx f 

Jo Jo 



dy- 



[( y + l-2x) 2 +y^} 2 



(65) 
(66) 



Appendix C 

Consider the contact potential given at any chiral order. In the matrix form defined in Sec. III. A, we have 
V = U T XU,T = U t tU, with E/(p) = (l,p 2 ,p 4 ,p 6 , . . .) being a column vector and U T being the transposed 
vector. Then the convolution in LSE could be factorized as VGqT — U 1 \p)XEtV \q) , with the matrix I 
being defined as follows, 



kdk 2 U{k)U T {k) 
(2tt) 2 E- k 2 /M + ie' 



(67) 
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The 3x3 case of I is given in Appendix A. It is easy to see that we could rewrite X as follows, 

J = I U(VMe)U t (VmE) + l([J m ];ME),mJ= 0, (68) 

where Iq and J m with ra^O are defined in Appendix A. Here I is a real matrix independent of Iq. From 
Eq.® it follows that, 

VG T = U T X{I UU T + X)tU = I VT + U t XJtU (69) 
Then, using the parametrization in Q, we find that, for on-shell momentum: 

t- 1 = y- 1 -g = y- 1 -i -g, (70) 



with g = u = g — Iq- Now comparing this with the following representation of T derived in Ref.|ll| 

using the relation between the on-shell T-matrix and on-shell JC-matrix, T _1 = K~ x + j^ip, we could find 
that, 



g = y- 1 - k- 1 + j , (71) 

that is, g must be a real number. But this real quantity is constructed with a complex T that contains the 
infinite iterations of the complex number I as given in l|7(J|l . That means Iq must cancel out in the infinite 
iteration, and hence must disappear in the real quantity g. This in turn implies that, J , as the real part of 
7o, does not appear in g. Finally these facts will lead to following from of T-matrix constructed with local 
potential: 

T -i _ Emc..},[J..})p 2t T £/Vi([C-UJ-])p 2i m (79 . 

- E^i([C..], [J...])pM 7 ° ~ EDjdC..], [J..])pV +J ° + 4 / P ' 
with [Ni,Dj] being independent of Jq. QED. 
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Fig.l Tree vertex for 4-nucleon in UT (left) and in EFT (right). 




Fig.2 Box diagram in UT (left) and bubble diagram in EFT (right). 
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